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Conducting steady-states of doped bilayer graphene have a non-zero sublattice pseudospin polar- 
ization. Electron-electron interactions renormalize this polarization even at zero temperature, when 
the phase space for electron-electron scattering vanishes. We show that because of the strength 
of interlayer tunneling, electron-electron interactions nevertheless have a negligible influence on the 
conductivity which vanishes as the carrier number density goes to zero. The influence of interactions 
is qualitatively weaker than in the comparable cases of single-layer graphene or topological insula- 
tors, because the momentum-space layer pseudo spin vorticity is 2 rather than 1. Our study relies 
on the quantum Liouville equation in the first Born approximation with respect to the scattering 
potential, with electron-electron interactions taken into account self-consistently in the Hartree-Fock 
approximation and screening in the random phase approximation. Within this framework the result 
we obtain is exact. 



I. INTRODUCTION 

The unique properties of graphene, the first truly two- 
dimensional material, have spurred a flood of research on 
fundamental physics and on technological possibilities^^ 
Monolayer graphene (MLG) has a non-Bravais honey- 
comb lattice structure with two triangular sublattices. 
The physical consequences of its linear 7r-band cross- 
ing at the Fermi level, described at low energies by a 
chiral massless Dirac k ■ p Hamiltonian, have been dis- 
cussed at length. Bilayer graphene (BLG) consists of 
two Bernal-stacked coupled graphene monolayers, result- 
ing in four inequivalent sites and four corresponding tt- 
bands. Its carriers are chiral but massive and charac- 
terized at low energies by a gapless parabolic spectrum 
displaying sublattice pseudospin-momentum locking. In 
BLG a band gap may be induced by a top gate^^l or hy 
dual gates which can vary the carrier density in the two 
layers independently.^ The theory of electronic trans- 
port in this unique and tunable 7r-band system has been 
investigated extensively!^^ Experimental studies have 
focused on transport EHEH magnetotransportj^^l a nd 
optics IHiHSSl Recent efforts that have succ eeded in manu- 
facturing and manipulating quantum dots p3E5] have been 
motivated by potential applications in quantum comput- 
ing. 

An extensive body of research has been devoted to 
electron-electron interactions in BLG.^^U Interactions 
in all forms of graphene are expected to be strong 
when the two-dimensional material is surrounded by low- 
re dielectrics. Progress in studying interaction effects 
has been improved by achieving samples with higher 
mobility,^ by studying suspended BLG samples that are 
not influenced by a substrate^ and by break throu ghs 
in fabricating samples with top and back gates! 95 ! 96 ! j n _ 
teractions are expected to become more impo rtant as 
one approaches the charge neutrality point In equi- 
libriu m int eractions in BLG lead to competin g g round 
states including a host of exotic states.^!!!! Theo- 
retical studies of electron-electron interactions in equilib- 



rium BLG^H23 have demonstrated, among other proper- 
ties, that screening and Friedel oscillations have different 
functional forms from MLG and 2DEGsPH5Sl 

The role of electron-electron interactions out of equilib- 
rium has not yet received attention. Given the wealth of 
research on transport, it is timely to address the influence 
of electron-electron interactions on the charge conductiv- 
ity of BLG. As in single-layer graphene, charge currents 
in BLG lead to a net pseudospin polarization. Interac- 
tions are therefore expected to renormalize the charge 
current and with it the pseudospin polarization. The 
question naturally arises of whether this polarization may 
be enhanced by interactions and produce observable ef- 
fects. It is important to understand whether the effect on 
the conductivity of non-equilibrium contributions to the 
interaction self-energy can be substantial, and whether it 
can be controlled using various tuning parameters such 
as the carrier density n e . 

This paper is therefore concerned with the effect of 
electron-electron interactions in bilayer graphene trans- 
port in the metallic regime £ft/H 3> 1. We begin with 
the quantum Liouville equation for the density matrix, 
working in the first Born approximation with respect 
to momentum scattering. Electron-electron interactions 
are taken into account self-consistently using the non- 
equilibrium Hartree-Fock approximation, with screening 
treated in the random phase approximation. We deter- 
mine an exact expression for the conductivity in the pres- 
ence of interactions within our framework. This work is 
distinct from recent papers discussing other interactions 
in BLG transport P3M1 j n addition, the mean-field ef- 
fect discussed here is not related to Coulomb drag, and 
electron-electron scattering is not relevant to the discus- 
sion at hand, which for simplicity assumes that the tem- 
perature T = 0. 

We demonstrate that electron-electron interactions 
renormalize the charge conductivity. The interaction ef- 
fect reduces the conductivity. However, the effect has 
a very weak density dependence and will be difficult to 
distinguish experimentally from a slight increase in disor- 
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der strength, which is not normally known. Surprisingly, 
the interaction effect vanishes as the carrier density n e 
tends towards zero because of a subtle interplay between 
the electric field, the pseudospin degree of freedom, and 
the electron-electron interactions mean-field. The effect 
is unexpectedly weak when the the Fermi wave vector hp 
is small compared to a wave vector introduced below, 
the size of which is set by the interlayer tunneling. In 
BLG this wave vector is go ~ 4nm _1 . Consequently even 
a density ~ 10 13 cm -2 , which is relatively large in experi- 
mental terms, gives a Fermi wave vector small compared 
to q a . 

A recent study of out-of-equilibrium interactions in TI 
showed that, for a Dirac cone, the renormalizcd con- 
ductivity has the same density dependence as the bare 
conductivity.^ This result is expected to apply also to 
MLG. In BLG, however, one would expect the renormal- 
ization to have a different density dependence given the 
quadratic dispersion and different functional form of the 
screened Coulomb potential. We find that the density de- 
pendence of the renormalization is indeed different from 
that of the bare conductivity, but the fractional change in 
the conductivity due to interactions is much weaker than 
in TI/MLG and vanishes at low densities. The vanishing 
of the renormalization in the limit n e — > is explained by 
the fact that the pseudospin of BLG is characterized by a 
winding number of 2. The projection of the equilibrium 
pseudospin at k onto the pseudospin at k' has a different 
rotational symmetry than the driving term due to the 
electric field. As n e — > 0, the product of these terms av- 
erages to zero over the Fermi surface. In this sense the 
relative weakness of this interaction effect in BLG is re- 
lated to the smaller Fermi velocity enhancement in BLG 
compared to MLGP 1 

The outline of this paper is as follows. In Sec. [TTJ we in- 
troduce the BLG band Hamiltonian and discuss the tech- 
nical details of the kinetic equation solution. In Sec. |III| 
we calculate the scattering term in BLG in the Born ap- 
proximation. In Sec. IV we briefly review charge trans- 



port in the absence of interactions. In Sec. [V] we calcu- 
late the first order mean-field correction to the conduc- 
tivity due to electron-electron interactions, then obtain 
an exact result to all orders. The results are discussed in 



Sec. VI while Sec. VII summarizes our findings. 



II. HAMILTONIAN AND METHOD 

Our study is based on a commonly used 2-band model 
for BLG: 



constant, which determines the Fermi velocity of BLG. 
Here v F is the (constant) Fermi velocity of MLG and 
tj_ w 0.4cV is the interlayer hopping parameter^ This 
model is valid at energies small compared to t±, except 
that it neglects trigonal warping terms which become im- 
portant at very low energies. (We comment on the role 
of these terms in Sec. VI ) It acts in a layer pseudospin 



space and has eigenstates which are equal weight sums 
of top and bottom layers with interlayer phase angle 29. 
The Hamiltonian Hg k can be understood as represent- 
ing a Zeeman-like interaction involving the pseudospin 
degree of freedom with a momentum-dependent effec- 
tive magnetic field whose direction is given by fi k . Un- 
like MLG, the pseudospin winds twice when momentum 
winds around the Fermi surface. The eigenvalues of H ok 
are e k ± — ±Ak 2 . 

The many-body Hamiltonian H in 2 nd quantization is 



H = 



E 



,t j 

1 L -k+q,s C k'-q,s' i 



'Cfe s ). 



(2) 

The one-particle matrix element H kk > ss i accounts for 
band structure contributions, as well as disorder and 
driving electric fields, discussed below. For the matrix 
element V q = V q we use the statically screened Coulomb 
potential, determined here in the random phase approx- 
imation (RPA)pSlwith e r the relative permittivity, 



Vn 



2e e r [q + q g(q)} ' 



(3) 



where the constant wave vector qg 



2ne a e r A 



and 



2k p 



44 



g 4 - In 



k 2 F + y/kjr + q 4 /4 



2k p 



(4) 



g(q) is a dimensionless function that increases monoton- 
ically from 1 to 1.755 as q varies from to 2k p, and the 
Fermi wave vector k F = y/Trn e . Both intra-band and 
inter-band contributions to static screening are included 
in g(q). For definiteness we will assume that the carrier 
density n e > 0. 

The effective single-particle kinetic equation for the k- 
diagonal part of the density matrix, f k , is derived from 
the quantum Liouville equation in the weak momentum 
scattering regime exactly as in Ref. |H7] 



Um = •!/••- ( e ie G q 6 I =Mr«T-n k .) (1) 



where fik is the unit vector (cos 20, sin26*) T with 9 the 
polar angle of k and A = h 2 v F /t± is a material-specific 



^ + \[H 0k ,f k ] + J(f k ) = -\[H*,f k ] + \[B^,f k ], 

(5) 

The scattering term J(f k ) in the first Born approxima- 
tion is given by 
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J(fk) = r| lim 



d 2 k' 



\U kk ,\ 2 / d^ e -^'{e- <ff ^^(/ fc _/ fc ,) e *ffo**7'' + e -i»o*t'/ft (/fc _ /fc/)ei H oh , t '/n }j (6) 



r 



where is impurity density and the potential of a 
single impurity. The mean-field electron-electron inter- 
action term is 



where Z = 1 is the ionic charge. Below we suppress the 
pseudospin indices ss' and treat all quantities as 2 x 2 
matrices. 



i 



(2tt)S 



dk'k' / e?7 V kk > fy 



(7) 



where 7 = 6' — 9 is the relative angle between wave vectors 
k and fc', Vfcfc' = V q , and and 9' are the polar angles 
of k and fc' respectively. The one-particle Hamiltonian 
HEkk'S SS ' + U kk >5 ss >, where H Ekk > is 



fefe' 



the electrostatic potential due to the driving electric field 
E, and U kk i is the foia^ disorder potential. The matrix 
element U kk > of the RPA-screened Coulomb potential of 
a single impurity between plane waves is 



U kk ' 



Ze 2 



2e e r [q + q a g(q)] 



(8) 



We decompose f k — n k TL + S k , with n k a scalar part 
and S k a pseudospin part, which can be expressed as 
S k = \S k ■ cr, where the vector S k is real and its z 
component is zero in equilibrium. Since the current 
operator is proportional to <x, we are only interested 
in S k . We decompose J(f k ) = J{n k ) + J(S k ) and 
{fk) = Bi IF (n k ) + F {S k ), and S k satisfies 



dS k 
dt 



[H ok , S k ] + J(S k ) = --[H£,S k ] + -[B™* (S k ),S k ] 



(9) 



The interaction terms in Eq. ([9} can be included itcra- 
tively, i.e. the solution is expanded in orders of V, i.e. 



S k - E„ S [ 



ee,(n) 
Ek 



MF 



0, and B 



MF,(r. 



B 



(2tt): 



V B 



MF,(n) 



with B 



MF,(Q) 



/ dk'K ft* dnv kk .s_ 



Ek 



,(n-l) 



Here 



5^fc = in the absence of electron-electron interac- 



tions. Substituting the above expansion into the kinetic 
equation and keeping only terms linear in the electric field 
(note that B^f F is linear in E), we obtain the equations 
below for each order n> Q: 



dS' 



ee,(n) 
Ek 

dt 



^[ H ok,S Ek 



(n)l 



J[S 



Ek 



[H k ,S[ 



Ek 



,(n)i 



(10) 



where Sofc is the pscudospin-dcpcndcnt part of the equi- 
librium density matrix, given below. 

Following the method used in Ref. we solve Eq. 



( 10 1 for each n by projecting onto directions parallel to 



(commuting with) and perpendicular to H$ kl obtaining 



where the parallel and perpendicular components 
S H = (l/2)(S fe • h k )(cT ■ h k ) = (l/2)s fc ||<T fc || 



S k _L = (1/2) (S fc • rh k )(cr ■ rh k ) 
and the unit vector rh k = z x n k . 



(l/2)sfe_LO- feJ _ 



(12) 



dt 



^ + ^[H k ,S kl ]+P ± J(S k ) 



P {] J(S k )= D H 
D k ±, 



III. SCATTERING TERM 

(H) The scattering term does not mix the monopole (n k ) 

and dipole (S k ) components of the density matrix. Ap- 
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plying the decomposition S k = Sfe||^fe + s k ±rh k and tr = 
&k\\'f l k + o'k±'riT'k, as well as rifc' = cos(27)nfc+sin(27)mfc 
and i7i/~> = — sin(27)nfc + cos(27)rh;,, to the expression 
for J(S k ), we obtain four projected terms as 



assume the temperature to be absolute zero, thus 



P\\J(S k \\) 
PxJ(S k \\) 
P\\J(S kX ) 
P±J(S k ±_) 



16tt/l4 
l6whA 
16nhA 
WttHA 



de'\Ukk'\ 2 {s k \\ -s fe '||)(l + cos27) 

de'\U kk ,\ 2 {s H - s fci )sin27 

d9'\Uhk'\ 2 (s k x +Sfe'x)sin27 

de'\U kk ,\ 2 {s kl _ + s fe / X )(l - cos 2 7 ) 

(13) 



Using q = 2kp sin 2r, we obtain a cumbersome expres- 
sion for I U kk > I . We make the following Fourier expansions 

\Ukk>\ 2 {i) = Y. u ™ eini 

(l + cos2 7 )|C7 fefe ,| 2 (7)=^W ne ^ (14) 

Sfc ll = s k\\n elnB ■ 

The parallel projection of the scattering term 

P\\J(S k \\) = ^^(Wo-W^a^e™*^ ( 15 ) 

n 

where W- n — W n since \ U kk i \ 2 {"f) and 1+cos 27 are even 
functions of 7. 



IV. TRANSPORT IN NON-INTERACTING BLG 

We briefly review transport in the absence of interac- 
tions. Writing S k — So k + Seic and keeping terms to o 
linear order in E we obtain 



— 37^ + T[Hok, Spk] + J(^Ek) — P>Eh 
at n 



(16) 



where the electric-field driving term 



D 



Ek 



l Ek\\ 



lEkX 



eE dSok 1 



dEk\\ a k\\ 



1 



h dk 2 

eE-k (df 0+ df - 
h \ dk dk 

eE - § (f f 1 
■l/o+ - Jo- J, 



■ dEk±Pk± 



(17) 



in which /o± = jb(£fc±), with /o the Fermi-Dirac dis- 
tribution function, and So k = (l/2)(/o+ — /o-)°fc||- Wc 



l Ek\\ 



5i. 



eE ■ k 

reE 



5(k 
k 



AH 



5(k - k F )a k \\ 



where the momentum relaxation time 



8hA 



riiiWo-Wi) 



The velocity operator is given by 



Vk = 



lMok 
h dk ' 



(18) 



(19) 



(20) 



The expectation value of the current density opera- 

/d 2 k 
j—yz Tv[v k Sk], where Tr acts 

in pseudospin space, and g v = g s = 2 are the val- 
ley and spin degeneracies, respectively. Substituting 
S k = (l/2)(s fe y(T fe || + s k ±a k ±), taking E || x, the 
velocity operator is v x = WzfeycrfcH + v xk ±a k ±, where 
v xk ii = 2Akcos9/h and v xk ± = —2Aksm6/h, and as- 
suming that £pT jh 1, it follows that the conductivity 
is 



(21) 



irk 2 



The Zitterbewegung (interband coherence) contribution 
to the conductivity, plays an essential role for the mini- 
mum conductiv ity which surviv es at t he charge neutrality 
pointj 15 | 2ti | 98 | 22l, as in_thLe MLGpEHfiland topological in- 
sulators (TI) cases 1121U22I but is next-to-leading order in 
the small parameter h/e F r and not considered here. 



V. INTERACTION RENORMALIZATION 

Interactions in equilibrium BLG renormalize the con- 
stant A (that is, they renormalize the Fermi velocity).^ 
This does not make any qualitative changes to our argu- 
ments and derivation below, and for simplicity we assume 
henceforth that A represents the renormalized A. In this 
section we will determine the_mean-ficld interaction cor- 



rection B^ F '^\ From Eq. ilOh it is evident that only 



the part of B^ F '^ cc a k ± contributes to the dynamics. 
We abbreviate I = k/k F and 



B 



MF,(1) 
k 



re 3 E r 



167reoe r fi. 
where the dimensionless quantities 



/W(Z,n e )sin 9<r k ±, 



(22) 



5 



2tt 



eh/ 



lixn e sin 7 sin(27) 



2tt 0rn e (l 2 + V 2 - 21V cos 7) + q g(l, V, 7) 



ff(M',7) = (i 2 + ^' 2 - 2ZZ' cos 7 ) 2 - In - + i >/4 + (Z 2 + Z' 2 - 211' cos 7 ) 2 



(23) 



r 



The driving term arising from contributes only 

to Sk±, and we easily find that 



,(i) _ TeE x q (1 



J Ekl. 



lW{l,n e )f smOa 



(24) 



I6hk 2 

An additional correction arises from the equation 

PuJ[S?$]=-P\\J[S%& ) ], (25) 

Taking this into account, the first-order correction to the 
diagonal conductivity is 

<£? = ^^ ^{neWlT (26) 

with/3(n e ) = (U 1 -U 3 )/{2U + 2U 2 -SU 1 -U 3 ), in which 
U n is the n-th Fourier coefficient of |E/fcfc'| 2 as defined in 

Eq. (14 1, and lil\n e ) = J* dllil\l,n e ). Notice that f3 
vanishes for momentum-independent (short-range) inter- 
actions. 



The angular structure of Iee(l,n e ) in Eq. (23) can be 



understood by noting that the electric field driving term 
is responsible for the factor of sin 7, while the factor of 
sin 27 arises from the projection of the pseudospin com- 
ponent parallel to k onto the pseudospin component par- 
allel to h! . For the masslcss Dirac cones of TI and MLG, 



where the (pseudo)spin winds around the Fermi surface 
only once, these terms (i.e. the electric-field driving term 
and the pseudospin projection) have the same rotational 
symmetry and reinforce each other. In BLG, the fact that 
the pseudospin winds twice around the Fermi surface is 
crucial, and makes the angular structure of this term en- 
tirely different from MLG and TI. As n e -t 0, lH\l,n e ) 
averages to zero over the Fermi surface. Its effect at small 
n e is therefore correspondingly small. In this context, it 
must also be noted that q is set by t±, the (sizable) in- 
terlayer hopping parameter, and that 50 > fcj? even at 
n = 10 13 cm~ 2 , which in transport ordinarily constitutes 
a large carrier density {hp s» 5.5 x 10 8 m _1 ). 

We retain only terms of linear order in the external 
electric field. Under these conditions, the following two 
equations are sufficient to obtain all higher order terms 
(n> 1), 



dS' 



ee,(n) 
Ek± 



dt 



In the higher 

,MF,(n) 



1 \vr c ee >(™)l 



MF,(r 



, Sok] 



(27) 



P \\J[S 
orders 



Ek 

(n 



> 1), s 



'Ek 



qee,(: 
^Ek 

ee,(n) 



is fed into 



completing the 



B k , which then determines S Ek 
self-consistent loop. Repeating the iteration, we obtain 
a general formula for j|™ (n e ) for n > 1, i.e. 



J&>(n e ) = (-^rrTj 



n 



n— 1 ~x 



n 



dli 

U 



2tt 



2ir 



dl r , 



2tv 



27T 



[cos 7, cos(27i) + /3(n e ) sin 7, sin^)] 



i=l 



^Jim e (l 2 + l 2 i+l - 2l t l i+ i cos7i) + qog(k,li+i,ji 

sin7»sin(27 ra ) 

^/im e {]? n + 1 - 2l n C0S7„) + q g(ln, 1, 7n) ' 
I 



(28) 



and the nth-order interaction correction to the conduc- Finally, the exact conductivity is 
tivity is a™ = [1 + P(n e )}{q /Vl^~e) n lie\n e )a^ rc . 



^barc 



1 + [! + £«)] 



y qo 



(29) 
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FIG. 1: Fractional change in the conductivity <J xx /a x ' x re 
(black) , and the parameter f3 (blue) , as a function of the car- 
rier density n e . 



We refer to a xx as the full conductivity, to distinguish 
it from the bare conductivity cr^ re . The appearance of 
f3(n e ) <g; 1 is related to the factor of 2-f appearing in the 
mean-field interaction term. 



VI. DISCUSSION 

In equilibrium electron-electron interactions renormal- 
ize the band parameter A.^Here we have obtained an ex- 
act result, within a self-consistent Hartree-Fock approx- 



imation, for the influence of interactions on the conduc- 
tivity of doped metallic BLG. Below we comment on the 
sign of the interaction renormalization, its size, and its 
density dependence. 

We recall that a charge current in BLG necessarily 
gives rise to a steady-state pseudospin polarization. Con- 
sequently, a xx may be understood by considering pseu- 
dospin dynamics on the Fermi surface. The renormaliza- 
tion reflects the interplay of pseudospin-momcntum lock- 
ing embodied in Hq^ and the mean pseudospin-field Bj^ IF 
arising from electron-electron interactions. The pseu- 
dospin of one carrier on the Fermi surface at k is subject 
to two competing interactions. The effective field Ak 2 hk 
tends to align the spin with its band value. The mean- 
field B^ IF tends to align a pseudospin at k against the 
total existing pseudospin polarization. The net result is 
a small steady-state rotation of the pseudospin at each 
k away from the direction of the effective field Ak 2 fif-. 
Thus the overall effect of interactions is to align individ- 
ual pseudospins in the direction opposite to that of the 
existing pseudospin polarization. 

Since the renormalization is negative, interactions can- 
not cause a xx to diverge, and there is no possibility of 
a Fermi-surface instability. The conductivity is therefore 
reduced by interactions, which is reminiscent of the result 
of Ref. 11031 One may gain insight by further analyzing 
the functional form of the ratio a xx /(j b x a x re 1 concentrating 
on its density dependence. Taking A = 0.71eV-nm 2 , 
as well as e r = 1 for simplicity, the wave vector qo = 
e 2 / '(2-Ke e r A) = 4.0nm -1 . As discussed before, in all 
realistic transport regimes kp — y/nn e <C qo- In this 
low-doping regime, /3(n e ) becomes independent of n e for 

1/2 

all U n oc n e , and the conductivity simplifies to 



T barc 



! + (!+/?)£ 



n 



1 dh 



2tt 



2^ 



dl, 



djn 

„ 2^ 



n 

,i=l 



cosji cos(27i) + /3sin7j sin(27i) 



sin7„sin(27 ra ) 
g(l n , 1,7„) 



(30) 



r 



In this limit the full conductivity has almost exactly the 
same density dependence as the bare conductivity. The 
behavior at densities commonly encountered in transport 
is illustrated in Fig. [l] The small size of the renormal- 
ization makes its detection challenging. At small n e , the 
ratio tends to zero as a result of the va nish ing of angular 



integral appearing in Iee(l,n e ) in Eq. (23). Steady-state 



expectation values are determined by the electric-field 



driving term, which contains a factor of E ■ k. Unlike 
MLG/TI, the pseudospin is not a linear function of || k 
(or 6), but is characterized by a winding number of 2. As 
a result of this, in BLG the interaction renormalization 
of the conductivity /pseudospin polarization is negligible 



when kp <C qo- At large n e 
marized by 



the behavior of a xx is sum- 
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(l+fl)go 
-y/167rn e 



1 /-27T 

dl 
o Jo 



g?7 sin7sin(27) 
2tt ^2 + 1-2/0037' 



(31) 



r 



In this regime the ratio cr xx /cr^!' e ^ VV^e' decreases 
with increasing carrier density, but for this trend to be 
noticeable one requires y/ivn e 3> (70 > which can never be 
reached in practice. 

It is enlightening to compare the interaction renormal- 
ization of the conductivity in BLG with the case of TI 
and MLG. In TI, as in MLG, the interaction renormal- 
ization of the conductivity is density independent and 
again accounts for only a fraction of the total conduc- 
tivity. At first sight, it seems striking that the same 
observation holds in BLG. Retracing the mathematical 
steps, the first order correction to the density matrix in 
TlisP 



,(i) _ eE x r s Tljl\l,r s ) 



S' 



EkJ_ 



/ sin 9(7 kJ _, 



(32) 



hence 1/k in TI corresponds to kpjk 2 in BLG, which 
results in approximately the same density dependence. 
The reason for this correspondence is that the TI Hamil- 
tonian is oc k while the BLG Hamiltonian is oc k 2 , so the 
steady-state (pseudo)spin densities differ by a factor of 
k. At the same time, screening also differs by a factor 
of k between the two, and the additional density depen- 
dences arising from these two factors effectively cancel 
out. Although the density dependence is different from 
TI, the correction is more complex but still weak. At 
very low-energies trigonal warping terms must be added 
to the BLG band structure, leading to the formation of 



Dirac cones. In this limit, we would expect that the inter- 
action correction to conductivity we discuss, would cross 
over to a form similar to that appropriate for MLG, TPs 
and other Dirac cone systems provided that this regime 
is not preempte d by in teraction-driven phase transitions 
to gapped statesPEsa 



VII. SUMMARY 

We have calculated the effect of non-equilibrium inter- 
action self-energy effects on the conductivity of metallic 
bilayer graphene. Although these effects can be large in 
some systems, in BLG they give rise to a negative renor- 
malization of the conductivity which is small and has a 
weak density dependence. This property follows from the 
large interlayer tunneling parameter in BLG, which leads 
to a 7r-band pseudo spin with a momentum-space wind- 
ing number of 2 that is incommensurate with the velocity 
winding number of 1. The corresponding effects could be 
larger when a gap is opened using a bias voltage or when 
a magnetic field is present. 
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